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Abstract 

This paper establishes the continuity of the path delay operators for dynamic network loading (DNL) problems based 
on the Lighthill-Whitham-Richards model, which explicitly capture vehicle spillback. The DNL describes and pre¬ 
dicts the spatial-temporal evolution of traffic flow and congestion on a network that is consistent with established 
route and departure time choices of travelers. The LWR-based DNL model is first formulated as a system of partial 
differential algebraic equations (PDAEs). We then investigate the continuous dependence of merge and diverge junc¬ 
tion models with respect to their initial/boundary conditions, which leads to the continuity of the path delay operator 
through the wave-front tracking methodology and the generalized tangent vector technique. As part of our analysis 
leading up to the main continuity result, we also provide an estimation of the minimum network supply without resort 
to any numerical computation. In particular, it is shown that gridlock can never occur in a finite time horizon in the 
DNL model. 
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1. Introduction 


Dynamic traffic assignment (DTA) is the descriptive modeling of time-varying flows on traffic networks consistent 
with traffic flow theory and travel choice principles. DTA models describe and predict departure rates, departure 
times and route choices of travelers over a given time horizon. It seeks to describe the dynamic evolution of traffic 
in networks in a fashion consistent with the fundamental notions of traffic flow and travel demand; see [Peeta and] 
ZiliaskopoulosI pOO 1 ) for some review on DTA models and recent developments. Dynamic user equilibrium (DUE) 
of the open-loop type, which is one type of DTA, remains a major modern perspective on traffic modeling that enjoys 
wide scholarly support. It captures two aspects of driving behavior quite well: departure time choice and route choice 
(Friesz et al. 1993\ . Within the DUE model, travel cost for the same trip purpose is identical for all utilized route- 
and-departure-time choices. The relevant notion of travel cost is effective travel delay, which is a weighted sum of 
actual travel time and arrival penalties. 

In the last two decades there have been many efforts to develop a theoretically sound formulation of dynamic 
network user equilibrium that is also a canonical form acceptable to scholars and practitioners alike. There are two 
essential components within the DUE models: (1) the mathematical expression of Nash-like equilibrium conditions; 


* Corresponding author 

Email addresses: k. han@imperial .ac.uk (Ke Han ), piccoliOcamden. rutgers. edu (Benedetto Piccoli), tf riesz@psu. edu (Terry L. 
Friesz) 


Preprint submitted to Transportation Research Part B 


March 22, 2016 






















ARTICLE LINK: http://www.sciencedirect.com/science/article/pii/S0191261515002039 
PLEASE CITE THIS ARTICLE AS 

Han, K., Piccoli, B., Friesz, T.L., 2015. Continuity of the path delay operator for dynamic network loading with 
spillback. Transportation Research Part B, DOI: 10.1016/j.trb.2015.09.009 


and (2) a network performance model, which is, in effect, an embedded dynamic network loading (DNL) problem. 
The DNL model captures the relationships among link entry flow, link exit flow, link delay and path delay for any 
given set of path departure rates. The DNL gives rise to the notion of path delay operator, which is viewed as a 
mapping from the set of feasible path departure rates to the set of path travel times or, more generally, path travel 
costs. 

Properties of the path delay operator are of fundamental importance to DUE models. In particular, continuity of 
the delay operators plays a key role in the existence and computation of DUE models. The existence of DUEs is 
typically established by converting the problem to an equivalent mathematical form and applying some version of 
Brouwer’s fixed-point existence theorem; examples include Han et al. (201 3c[>; Smith and Wistenjl 1 995[>; |Wie et al~ 


( 2002) and Zhu and Marcotte ( 2000| > . All of these existence theories rely on the continuity of the path delay operator. 
On the computational side of analytical DUE models, every established algorithm requires the continuity of the delay 
operator to ensure convergence; an incomplete list of such algorithms include the fixed-point algorithm (Friesz et al., 
|2013| >, the route-swapping algorithm ( Huang and Lam|[2002| l, the descent method ( |Han and Lo| 2003), the projection 
method (Han and Lo 2002| Ukkusuri et al. 2012| >, and the proximal point method (Han et al. 2015a| > 

It has been difficult historically to show continuity of the delay operator for general network topologies and traffic 
flow models. Over the past decade, only a few continuity results were established for some specific traffic flow 
models. Zhu and Marcotte ( 2000) > use the fink delay model (Friesz et al. 19931 to show the continuity of the path 
delay operator. Their result relies on the a priori boundedness of the path departure rates, and is later improved by 
a continuity result that is free of such an assumption ( |Han et aT] |2012| >. In|Bressan and Han|(|2013|l, continuity of 


the delay operator is shown for networks whose dynamics are described by the LWR-Lax model (Bressan and Han 


2011 [ Friesz et al. 20131, which is a variation of the LWR model that does not capture vehicle spillback. Their result 


also relies on the a priori boundedness of path departure rates. Han et al. ( |2013c I consider Vickrey’s point queue 
model ( Vickrey) 19691 and show the continuity of the corresponding path delay operator for general networks without 
invoking the boundedness on the path departure rates. 

All of these aforementioned results are established for network performance models that do not capture vehicle 
spillback. To the best of our knowledge, there has not been any rigorous proof of the continuity result for DNL 
models that allow queue spillback to be explicitly captured. On the contrary, some existing studies even show that 
the path travel times may depend discontinuously on the path departure rates, when physical queue models are used. 
For example, Szeto ( |2003| > uses the cell transmission model and signal control to show that the path travel time may 
depend on the path departure rates in a discontinuous way. Such a finding suggests that the continuity of the delay 
operator could very well fail when spillback is present. This has been the major hurdle in showing the continuity or 
identifying relevant conditions under which the continuity is guaranteed. This paper bridges this gap by articulating 
these conditions and providing accompanying proof of continuity. 

This paper presents, for the first time, a rigorous continuity result for the path delay operator based on the LWR 
network model, which explicitly captures physical queues and vehicle spillback. In showing the desired continuity, 


a diverge model, both originally presented by 

)aganzo (1995). The underpinning analytical framework employs 

the wave-front tracking methodology (Dafermos 

1972 Holden and Risebro 2002) and the technique of generalized 

tangent vectors (iBressanj 1993; Bressan et al. 

2000). A major portion of our proof involves the analysis of the 


interactions between kinematic waves and the junctions, which is frequently invoked for the study of well-posedness 
of junction models; see Garavello and Piccoli ( 2006j ) for more details. Such analysis is further complicated by the 
fact that vehicle turning ratios at a diverge junction are determined endogenously by drivers’ route choices within the 
DNL procedure. As a result, special tools are developed in this paper to handle this unique situation. 

As we shall later see, a crucial step of the process above is to estimate and bound from below the minimum network 
supply, which is defined in terms of local vehicle densities in the same way as in Lebacque and Khoshyaran ( j 1 999) . 
In fact, if the supply somewhere tends to zero (that is, when traffic approaches the jam density), the well-posedness of 
the diverge junction may fail, as we demonstrate in Section |4.2.1 This has also been confirmed by the earlier study 
of Szeto ( 2003| >, where a wave of jam density is triggered by a signal red light and causes spillback at the upstream 
junction, leading to a jump in the path travel times. Remarkably, in this paper we are able to show that (1) if the supply 


'Well-posedness of a model refers to the property that the behavior of the solution hardly changes when there is a slight change in the ini- 
tial/boundary conditions. 
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is bounded away from zero (that is, traffic is bounded away from the jam density), then the diverge junction model is 
well posed; and (2) the desired boundedness of the supply is a natural consequence of the dynamic network loading 
procedure that involves only the merge and diverge junction models we study here. This is a highly non-trivial result 
because it not only plays a role in the continuity proof, but also implies that gridlock can never occur in the network 
loading procedure in any finite time horizon. 

The final continuity result is presented in Section |5~4| following a number of preliminary results set out in previous 
sections. Although our continuity result is established only for networks consisting of simple merge and diverge nodes, 
it can be extended to networks with more complex topologies using the procedure of decomposing any junction into 
several simple merge and diverge nodes (Daganzo 1995] ). Moreover, the analytical framework employed by this 
paper can be invoked to treat other and more general junction topologies and/or merging and diverging rules, and the 
techniques employed to analyze wave interactions will remain valid. 

The main contributions made in this paper include: 


• formulation of the LWR-based dynamic network loading (DNL) model with spillback as a system of partial 
differential algebraic equations (PDAEs); 

• a continuity result for the path delay operator based on the aforementioned DNL model; 

• a novel method for estimating the network supply, which shows that gridlock can never occur within a finite 
time horizon. 


The rest of this paper is organized as follows. Section[2]recaps some essential knowledge and notions regarding 
the LWR network model and the DNL procedure. Section p]articulates the mathematical contents of the DNL model 
by formulating it as a PDAE system. Section [4] introduces the merge and diverge models and establishes their well- 
posedness. Section^] provides a final proof of continuity and some discussions. Finally, we offer some concluding 
remarks in Section[6J 


2. LWR-based dynamic network loading 

2.1. Delay operator and dynamic network loading 

Throughout this paper, the time interval of analysis is a single commuting period expressed as [0, T] for some 
T > 0. We let P be the set of all paths employed by travelers. For each path p G V we define the path departure rate 
which is a function of departure time t 6 [0, T]\ 


h p (■) : [0, T] -> R + 

where R+ denotes the set of non-negative real numbers. Each path departure rate h p (t ) is interpreted as a time-varying 
path flow measured at the entrance of the first arc of the relevant path, and the unit for h p (t) is vehicles per unit time. 
We next define h(-) = {h p (•) : p G V) to be a vector of departure rates. Therefore, h = /;(•) can be viewed as a 
vector-valued function of t, the departure time^j 

The single most crucial ingredient t is the path delay operator, which maps a given vector of departure rates h to a 
vector of path travel times. More specifically, we let 

D p (t, h) it G [f 0 , t f ], ipeP 

be the path travel time of a driver departing at time t and following path p, given the departure rates h associated with 
all the paths in the network. We then define the path delay operator D(-) by letting Dili) = \l) p (-, h) : p G P}, which 
is a vector consisting of time-dependent path travel times D p (t , h). 
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For notation convenience and without causing any confusion, we will sometimes use h instead of h{-) to denote path flow vectors. 
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2.2. The Lighthill-Whitham-Richards model on networks 

We recap the network extension of the LWR model (Lighthill and Whitham 1955; Richards] |1956) , which captures 
the formation, propagation, and dissipation of spatial queues and vehicle spillback. Discussion provided below relies 
on general assumptions on the fundamental diagram and the junction model, and involves no ad hoc treatment of flow 
propagation, flow conservation, link delay, or other critical model features. 

We consider a road link expressed as a spatial interval [a, b] c R. The partial differential equation (PDE) 
representation of the LWR model is the following scalar conservation law 

d r p(t, x) + d x f{p{t , x)) = 0 (f, x) e [0, T] x [a, b] (2.1) 


with appropriate initial and boundary conditions, which will be discussed in detail later. Here, p(t, x) denotes vehicle 
density at a given point in the space-time domain. The fundamental diagram /(•) : [0, p-' 1 ""] —> [0, C] expresses 
vehicle flow at (f, x) as a function of p(t, x), where p /V "" denotes the jam density, and C denotes the flow capacity. 
Throughout this paper, we impose the following mild assumption on /(•): 

(F). The fundamental diagram /(•) is continuous, concave and vanishes at p — 0 and p = p Jum . 


An essential component of the network extension of the LWR model is the specification of boundary conditions 
at a road junction. Derivation of the boundary conditions should not only obey physical realism, such as that enforced 
by entropy conditions (Garavello and Piccoli |2006 Holden and Risebro 1995[ >, but also reflect certain behavioral 
and operational considerations, such as vehicle turning preferences (Daganzo 1995), driving priorities (Coclite et 


ah] [2005), and signal controls (Han et aTj 201 4). Articulation of a junction model is facilitated by the notion of 


Riemann Problem , which is an initial value problem at the junction with constant initial condition on each incident 
link. There exist a number of junction models that yield different solutions of the same Riemann Problem. In one line 
of research, an entropy condition is defined based on a minimization/maximization problem (Holden and Risebro 
1995). In another line of research, the boundary conditions are defined using link demand and supply (Lebacque 


and Khoshyaran 1999), which represent the link’s sending and receiving capacities. Models following this approach 
include Daganzop 1995| ); [JiiTa nd Zhang (2003) and |Jin| ( f2010[ ). The solution of a Riemann Problem is given by the 
Riemann Solver (RS), to be defined below. 


2.2. 1. The Riemann Solver 

We consider a general road junction J with m incoming roads and n outgoing roads, as shown in Figure |T| 

m m+n 



Figure 1: A road junction with m incoming links and n outgoing links. 

We denote by the incoming links and by I m+ 1 ,..., I m+n the outgoing links. In addition, for every i e 

{1,... ,m + n), the dynamic on /, is governed by the LWR model 

d,pi(t, x) + d x fj(pi(t, x)) = 0 (f, x) e [0, T] x [a,-, b,] (2.2) 

where link Ij is expressed as the spatial interval [a,, /;,], and we always use the subscript T to indicate dependence on 
the link /,. The initial condition for this conservation law is 


Pi( 0, x) = pi(x) 
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X € [fl;, /?;] 


(2.3) 
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Notice that the above ( m+n) initial value problems are coupled together via the boundary conditions to be specified 
at the junction. Such a system of coupling equations is commonly analyzed using the Riemann Problem. 

Definition 2.1. (Riemann Problem) The Riemann Problem at the junction J is defined to be an initial value problem 
on a network consisting of the single junction J with m incoming links and n outgoing links, all extending to infinity, 
such that the initial densities are constants on each link: 


lp;(0, x) = pi x £ (- 00 , bj], i £ {1, ..., m\ 

lp ; (0, x ) = pj x £ [aj , + 00 ), j e [m +1, m + n\ 

where pk £ [0, p J “ n ] are constants, k — 1,... ,m + n. 

A Riemann Solver for the junction / is a mapping that, given any (m + n)-tuple of Riemann initial conditions 
(pi,... ,p m+ „), provides a unique (m + n)-tuple of boundary conditions (p 1; ... ,p m+ „) such that one can solve the 
initial-boundary value problem for each link, and the resulting solutions constitute a weak entropy solution of the 
Riemann Problem at the junction^] A precise definition of the Riemann Solver is given as follows. 

Definition 2.2. (Riemann Solver) A Riemann Solver for the junction J with m incoming links and n outgoing links is 
a mapping 


RS : Y] [ 0 , pf] p| [ 0 . pf"| 

k =1 k =1 

(Pl- • • ■ . Pm+n) l—> (Pl> • ■ • > Pm+n ) 

which relates Riemann initial data p = (pi,... ,p m+ „) to boundary conditions p = (p l ,..., p m+n ), such that the 
following hold. 

(i) The solution of the Riemann Problem restricted on each link f is given by the solution of the initial-boundary 
value problem with initial condition p k and boundary condition p k , k — 1,... ,m + n. 

(ii) The Rankine-Hugoniot condition (flow conservation) holds: 

m m+n 

Yjfifid - 2 (2.4) 

i— 1 j—m +1 

( Hi) The consistency condition holds: 

RS[RS\p]] = RSlp] (2.5) 


Three conditions must be satisfied by the Riemann Solver (RS). Item ( i ) above requires that the boundary condition 
on each link must be properly given so that the initial value problems not only have well-defined solutions, but these 
solutions must also be compatible and form a sensible solution at the junction. \2A) simply stipulates the conservation 
of flow across the junction. ( |2.5[ ) is a desirable property and is sometimes referred to as the invariance property (Jin 

2'OTO ). 

Remark 2.3. For the same Riemann Problem, there exist many Riemann Solvers that satisfy conditions (i)-(iii) above. 
Despite their varying forms, most existing Riemann Solvers rely on a flow maximization problem at the relevant 
junction subject to constraints related to turning ratio, right-of-way, or signal controls; see ( {Cocl ite et al. 
et al. \2014 Holden and Risebro 1995 ■ Jin and Zhang 2003) and Jin (2010). 


2005 • Han 


' W’e refer the reader to 


Holden and Risebro 


1995 for a definition of weak entropy solution at a junction 
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2.2.2. The link demand and supply 

For each link /„ we let // be the critical density at which the flow is maximized. The demand D,(l) for incoming 
links /, and the supply S jit) for outgoing links Ij are defined in terms of the density near the exit and entrance of the 
link, respectively (Lebacque and Khoshyaran 1999] ): 


D,(t) = D{ Pi (t, b,-)) 

Sj(t) = Sj(pj(t , aj+)) 


Ci 

-)) 

if P iit, b~) > p 1 : 

(2.6) 

fiipiit , bj- 

if piit, bi-) < p\ 

Cj 

''+)) 

if pjit, aj+) < p c j 

(2.7) 

fjiPjit’ 

if pjit, aj+) > p c j 


In prose, the demand represents the maximum flow at which cars can be discharged from the incoming link; and the 
supply represents the maximum flow at which cars can enter the outgoing link. Notice that the demand and supply 
are both expressed as functions of density, and they are always greater than or equal to the fundamental diagram 
f(-) or /,(•); see Figure [2] for an illustration. In our subsequent presentation, without causing confusion we will use 
notations D,(f) and D,(p) interchangeably where the former indicates the demand as a time-varying quantity, and the 
latter emphasizes demand as a function of density. The same applies to the supply. 



fip) 

Dip) 

Sip) 


3. Dynamic network loading problem formulated as a PDAE system 

The aim of this section is to formulate the LWR-based dynamic network loading (DNL) problem as a system of 
partial differential algebraic equations (PDAEs). The proposed PDAE system uses vehicle density and queues as the 
primary unknown variables, and computes link dynamics, flow propagation, and path delay for any given vector of 
departure rates. The PDAE system captures vehicle spillback explicitly, and accommodates a wide range of junction 
types and Riemann Solvers. 

We consider a network G(t7l, 'V) expressed as a directed graph with JK being the set of links and 'V being the set 
of nodes. Let V be the set of paths employed by travelers, and TV be the set of origin-destination pairs. Each path 
p eP is expressed as an ordered set of links it traverses: 

P — \J 1; Di - - ■ 11m(p)} 

where m(p) is the number of links in this path. There are several crucial components of a complete network loading 
model, each of which is elaborated in a subsection below. Throughout the rest of this paper, for each node v e T 7 , we 
denote by I v the set of incoming links and O' the set of outgoing links. 
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3.1. Within-link dynamics 

For each /, e J[. the density dynamic is governed by the scalar conservation law 

d,pi(t, x ) + d x |p,(f, x) ■ x))] = 0 (t, x) e [0, T] x [a h b,] 


(3.8) 


subject to initial condition and boundary conditions to be determined in Section 3.2 The fundamental diagram 
f(pi) - pi ■ \’i{pi) satisfies condition (F) stated at the beginning of Section 2.2 In order to explicitly incorporate 


drivers’ route choices, for every p e P such that /, e p we introduce the function p p (t, x), (t. x) e [0, T] x [a,-, Ip], 
which represents, in every unit of flow /)(p,(f, x)), the fraction associated with path p. We call these variables path 
disaggregation variables (PDV). For each car moving along the link /,, its surrounding traffic can be distinguished 
by path (e.g. 20% following path p \, 30% following path pi. 50% following path pi). As this car moves, such a 
composition will not change since its surrounding traffic all move at the same speed under the first-in-first-out (FIFO) 
principle (i.e. no overtaking is allowed). In mathematical terms, this means that the path disaggregation variables, 
p p {-, ■), are constants along the trajectories of cars (t, x(t)) in the space-time diagram, where x(-) is the trajectory of a 
moving car on link . That is, 

— u p (t, x(t)) = 0 ip such that /, e p, 
dt 

which, according to the chain rule, becomes 


x(t)) + djjJ-it, x ) • — x(t) = 0, 
dt 

which further leads to another set of partial differential equations on link /, : 

d t p p (t, x) + \’i(pi(t, x)) • d xt u p (t, x) = 0 ip such that /, 6 p 
Here, p,(f, x) is the solution of ( |3.8| |. The following obvious identity holds 

'y' p p (t, x) = 1 whenever p,(t, x) > 0 

P 3li 


(3.9) 


(3.10) 


where p 3 /, means “path p contains (or traverses) link and the summation appearing in ( |3. 1 ()[ > is with respect to 
all such p. By convention, if p,(f, x) = 0, then p p (t, x) = 0 for all p 3 /,. 

3.2. Boundary conditions at an ordinary node 

For reason that will become clear later, we introduce the concept of an ordinary node. An ordinary node is neither 
the origin nor the destination of any trip. We use the notation 'V 0 to represent the set of ordinary nodes in the network. 

As mentioned earlier, the partial differential equations on links incident to J are all coupled together through 
a given junction model, i.e., a Riemann Solver. A common prerequisite for applying the Riemann Solver is the 
determination of the flow distribution (turning ratio) matrix (Coclite et al. 2005j >, which relies on knowledge of the 
PDVs p p (t, bi) for all /, e I J . We define the time-dependent flow distribution matrix associated with J: 


A J (t ) = (ofTO) e [0, 1] |J 'I + I° J I 


(3.11) 


where by convention, we use subscript i to indicate incoming links, and j to indicate outgoing links. Each element 
ar/.(t) represents the turning ratios of cars discharged from /, that enter downstream link / ; . Then, for all p that traverses 
/, the following holds. 


a-j(t) = ^ p P (t, bd V I/ e I J 

pBli.Ij 


Ij 6 0 J 


(3.12) 


It can be easily verified that crjj(t) e [0, 1] and ; a/.(f) = 1 according to ( |3. 1 ()[ >. 

We are now ready to express the boundary conditions for the ordinary junction J e 'V" 


Let 


\l J \+\p J \ 


RS Al : f[ [0,pf m ] 


k= 1 


|J y l+l^| 

FI [°>pn 


k= 1 
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be a given Riemann Solver. Notice that the dependence of the Riemann Solver on A J has been indicated with a 
superscript. The boundary conditions for PDEs ( |3.8| ) read 

Pk(f, b k ) = RSf b,-)) Ii€jJ , (pj(t , a j+ )) Ij€OJ ] V4 e I J (3.13) 

pi(t, a t ) = RSf 1 bj-)), i€l j , cij+)) Ije0 j\ V/, e 0 J (3.14) 

where /?S^[-] denotes the k-th component of the mapping, k — 1, ..., \I J \ + \0 ] \. 

Remark 3.1. Intuitively, ( |3.13| )-( |3.14| i mean that, given the current traffic states (pft, bj—)) re jj and (pj(t, cij+)) I e0 i 

adjacent to the junction J, the Riemann Solver RS aJ specifies, for each incident link I k or the corresponding bound¬ 
ary conditions p k (t, b k ) or p/(t, at). In prose, at each time instance the Riemann Solver inspects the traffic conditions 
near the junction, and proposes the discharging (receiving) flows of its incoming (outgoing) links. Such a process 
is based on the flow distribution matrix A J and often reflects traffic control measures at junctions. Furthermore, the 
Riemann Solver operates with knowledge of every link incident to the junction, thus the boundary condition of any 
relevant link is determined jointly by all the links connected to the same junction. Therefore, the LWR equations on 
all the links are coupled together through this mechanism. For this reason the LWR-based DNL models are highly 
challenging, both theoretically and computationally. 


The upstream boundary conditions associated with PDEs (|3.91> are: 


P P j(t, af) 


fiipiit, bjj) • p p (t, bA 
fj(pj(t, aj )) 


'Ip such that {/,-, Ij) c p. 


Tlj e O' 


(3.15) 


where the numerator fiipft , /?,)) ■ p P (t, bj) expresses the exit flow on link /, associated with path p, which, by flow 
conservation, is equal to the entering flow of link Ij associated with the same path p\ the denominator represents the 
total entering flow of link I j. 

Remark 3.2. Unlike the density-based PDF, the PDV-based PDF does not have any downstream boundary condition 
due to the fact that the traveling speeds of the PDVs are the same as the car speeds (they can be interpreted as 
Lagrangian labels that travel with the cars); thus information regarding the PDVs does not propagate backwards or 
spills over to upstream links. 


3.3. Flow distribution at origin or destination nodes 

We consider a node v e 'V that is either the origin or the destination of some path p. One immediate observation 
is that the flow conservation constraint ( |2.4[ > no longer holds at such a node since vehicles either are ‘generated’ (if v 
is an origin) or ‘vanish’ (if v is a destination). A simple and effective way to circumvent this issue is to introduce a 
virtual link. A virtual link is an imaginary road with certain length and fundamental diagram, and serves as a buffer 
between an ordinary node and an origin/destination; see Figure [3] for an illustration. By introducing virtual links to 
the original network, we obtain an augmented network G(3\, 'V) in which all road junctions are ordinary, and hence 
fall within the scope of the previous section. 


Figure 3: Illustration of the virtual links. Left: a virtual link connecting an origin (s) to an ordinary junction J. Right: a virtual link connecting a 
destination (t) to an ordinary junction J. 




Let us denote by N the set of origins in the augmented network G(3\, 'V). For any s e S, we denote by V s c V the 
set of paths that originate from s, and by I s the virtual link incident to this origin. For each p e V s we denote by h p (t) 
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the departure rate (path flow) along p. It is expected that a buffer (point) queue may form at s in case the receiving 
capacity of the downstream I s is insufficient to accommodate all the departure rates f P eP’ h p (t). For this buffer queue, 
denoted q s (t ), we employ a Vickrey-type dynamic ([Vickrey). [1969}>; that is. 


^ h p (t) ■ 

peP’ 


S s (t ) 

min { ^ h p (t), 5 S (0} 


if q s (t) > 0 
if q s (t) = 0 


(3.16) 


pep’ 


where S s (t) denotes the supply of the virtual link I s . The only difference between < |3. 1 6[ ) and Vickrey’s model is the 
time-varying downstream receiving capacity provided by the virtual link. Q 

It remains to determine the dynamics for the path disaggregation variables (PDV). More precisely, we need to 
determine p p s {t, a s ) for the virtual link I s where p e V s , and x = a s is the upstream boundary of I s . This will be 
achieved using the Vickrey-type dynamic ( |3.16| ) and the FIFO principle. Specifically, we define the queue exit time 
function A s (f) where t denotes the time at which drivers depart and join the point queue, if any; Aft) expresses the 
time at which the same group of drivers exit the point queue. Clearly, FIFO dictates that 


r v- r AM 

V h p (r) dr = 

Jo ^ Jo 


f s (p s (T, a s )) dr 


(3.17) 


where the two integrands on the left and right hand sides of the equation are flow entering the queue and flow leaving 
the queue, respectively. We may determine the path disaggregation variables as: 

rf(Uf),a s ) = p(f) ip e P, (3.18) 

2uqep s ™q \0 

Notice that, if Yjqep 1 h q (t) = 0, then the flow leaving the point queue at time Aft) is also zero; thus there is no need to 
determine the path disaggregation variables. Therefore, the identity (|3. 1 8|) is well defined and meaningful. 


3.4. Calculation of path travel times 

With all preceding discussions, we may finally express the path travel times, which are the outputs of a complete 
DNL model. The path travel time consists of link travel times plus possible queuing time at the origin. Mathematically, 
the link exit time function Aft) for any /, is defined, in a way similar to (|3.17|i, as 




,(t, bi))dr 


(3.19) 


For a path expressed as p — {7), I 2 , ..., I m (p)}, the time to traverse it is calculated as 


A s o Ai o A 2 . ■. o A m (p)(t) (3.20) 

where / o git) = g{f{t)) means the composition of two functions. This is due to the assumption that cars leaving the 
previous link (or queue) immediately enter the next link without any delay. 

3.5. The PDAE system 

We are now ready to present a generic PDAE system for the dynamic network loading procedure. Let us begin by 
summarizing some key notations. 


4 The right hand side of the ordinary differential equation (3.16) is discontinuous. An analytical treatment of this irregular equation is provided 
by|Han et al.|(2013a|b) using the variational formulation. 
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Gift, °V) 

the original network with link set ft and node set 'V; 

A'Z 

the set of virtual links; 

Gift, V) 

the augmented network including virtual links; 

S 

the set of origins in G(ft, r V); 

V s 

the set of paths originating from s 6 S', 

<yo 

the set of ordinary junctions in Gift, 'V)\ 


the set of incoming links of a junction J e 'V "; 

o' 

the set of outgoing links of a junction J e 'V "; 

A J it) 

the flow distribution matrix associated with junction ./; 

rs aJ 

the Riemann Solver for junction J, which depends on A 1 . 


We also list some key variables of the PDAE system below. 


hpit) 

the path departure rate along p e V\ 

pAt, x) 

the vehicle density on link 7/ 6 ft. 

P P i it, x) 

the proportion of flow on link 7, associated with path p (path disaggregation variable); 

q s it) 

the point queue at the origin s e S', 

Ut) 

the point queue exit time function at origin .v e S. 


Given any vector of path departure rates h = : p e V), the proposed PDAE system for calculating path travel 

times is summarized as follows. 


dq s (t ) 
dt 


= X h P (t) ~ 


pef” 


S s (t ) q s (t)>0 

min iz* P (t), 5,(1)} qAA) = 0 


peP s 


f J] h p ( T ) dr = f fsiflsij, a s )) dT 

Jo Jo 

J rt 

fipiiT, dij) dr = fiipAr, bj) dT 

0 Jo 

dtPiit, x) + d x | pAt, x ) • ViipAt, a))] = 0 

x) + ViipAt, x)) ■ d x p p {t, x) = 0 

h p (t) 


p P MAt), a s ) 


P p j(t, aj) 


hqif) 

flip At, bj)) ■ p p (t, bi) 


fjipjit, aj )) 
A J (t) = [a]jit)}, aijit) = 2 p p (t, bj 


P3li,lj 


Pk(t, b k ) = RS f [ {pAt, bi-)) Ij£l j , {pj{t, «/+))/ (€ o'] 
pAt, a,) = RSf [ipAt, b~)) Ii€l j , ipjit, aj+)) Ije0 j] 
D p it, h) = A s o di o A 2 ... o A m ( P )(t) 


VseS 

VseS 
V/, e ^ 

it, x) e [0, T] x [a,, bj 
it, ,r) e [0, T] x [a,-, bi] 

Vi eS, p e Vs 

dp ^ {//- Iji 
V/, e I J , Ij e O' 

V4 e I J 
V/, e <9 7 

V p eV, dt e [0, T] 


(3.21) 

(3.22) 

(3.23) 

(3.24) 

(3.25) 

(3.26) 

(3.27) 

(3.28) 

(3.29) 

(3.30) 

(3.31) 
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Eqn p.21[ i describes the (potential) queuing process at each origin. Eqns ( | 3.22 [ ) and ( | 3.23 [ > express the queue exit 
time function for a point queue, and the link exit time function for a link, respectively. Eqns ( |3.24[ )-( [3T25| ) express the 
link dynamics in terms of car density and PDV; Eqns ( |3.26| l-( [T27| ) specifies the upstream boundary conditions for the 
PDV as these variables can only propagate forward in space. Eqns ( [3.28| >-( |3T3()| determine the boundary conditions at 
junctions. Finally, Eqn ( |3.31| > determines the path travel times. 

The above PDAE system involves partial differential operators d, and d x . Solving such a system requires solu¬ 
tion techniques from the theory of numerical partial differential equations (PDE) such as finite difference methods 
(Godunov 1959[ LeVeque 1992) and finite element methods (Larsson and Thomee 2005]l. 


4. Well-posedness of two junction models 

In mathematical modeling, the term well-posedness refers to the property of having a unique solution, and the 
behavior of that solution hardly changes when there is a slight change in the initial/boundary conditions. Examples 
of well-posed problems include the initial value problem for scalar conservation laws (Bressan 2000 [ >, and the initial 
value problem for the Hamilton-Jacobi equations (Daganzo 2006) >. In the context of traffic network modeling, well- 
posedness is a desirable property of network performance models capable of supporting analyses and computations 
of DTA models. It is also closely related to the continuity of the path delay operator, which is the main focus of this 
paper. 

This section investigates the well-posedness (i.e. continuous dependence on the initial/boundary conditions) of 
two specific junction models. These two junctions are depicted in Figure[4] and the corresponding merge and diverge 
rules are proposed initially by Daganzo ( 1995} in a discrete-time setting with fixed vehicle turning ratios and driving 
priority parameters. 




Figure 4: The diverging (left) and merge (right) junctions. 


4.1. The two junction models 

4.1.1. The diverge model 

We first consider the diverge junction shown on the left part of Figure [4] with one incoming link 7) and two 
outgoing links I 2 and It,. The demand D\(t) and the supplies, S 2 (f) and S 3 (f), are defined by ( | 2 .6 [ ) and ( |2.7[ ) respectively. 
The Riemann Solver for this junction relies on the following two conditions. 

(Al) Cars leaving 7) advance to It and 7 3 according to some turning ratio which is determined by the PDV p\ (t, h\) 
in the DNL model. 


(A2) Subject to (Al), the flow through the junction is maximized. 

In the original diverge model (Daganzo, 1995} , the vehicle turning ratios, denoted ct\T and 0-13 with obvious 
meaning of notations, are constants known a priori. This is not the case in a dynamic network loading model since 
they are determined endogenously by drivers’ route choices, as expressed mathematically by Eqn ( |3.28| >. The diverge 
junction model, described by (Al) and (A2), can be more explicitly written as: 


/r« = 

fi n (t) = 
f] n (t) - 


min 



Si(t) 

£*1,2(0’ 


*1,2(0-/no 

£*1,3(0-/no 


s 3 (0 ] 

£*1.3(0 J 
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where f° u, (t ) denotes the exit flow of link I \, and f"'(t) denotes the entering flow on link /,, j = 2, 3. 

4.1.2. The merge model 

We now turn to the merge junction in Figure [4] with two incoming links / 4 and I 5 and one outgoing link I h . In 
view of this merge junction, assumption (Al) becomes irrelevant as there is only one outgoing link; and assumption 
(A2) cannot determine a unique solution^] To address this issue, we consider a right-of-way parameter p e (0, 1) and 
the following priority rule: 

(Rl) The actual link exit flows satisfy (1 - p) ■ ff u, (t) = p ■ 



Figure 5: Illustration of the merge model. The shaded areas represent the feasible domain of the maximization problem j4.33) , and the thick line 
segments (£1) are the optimal solution sets of 14.33f Left: Rule (Rl) is compatible with (A2); and there exists a unique point Q satisfying both 
(A2) and (Rl). Right: (Rl) is incompatible with (A2); in this case, the model selects point Q within the set Cl that is closest to the ray (R) from the 
origin with slope —f. 


Notice that (Rl) may be incompatible with assumption (A2); and we refer the reader to Figure [5] for an illustration. 
Whenever there is a conflict between (Rl) and (A2), we always respect (A2) and relax (Rl) so that the solution is 
chosen to be the point that is closest to the line y = — x among all the points yielding the maximum flow. Clearly, 
such a point is unique. Mathematically, we let 11 be the set of points ( //'"', / 5 °"') that solves the following maximization 
problem: 

fma x/r + /r , 433) 

\ such that 0 < /“" < D 4 (t), 0 < f° ut < D 5 (t), f™ 1 + f° m < S 6 (t) y ’ 

Moreover, we define the ray R = {( ff l “, /“") e : (1 - p) ■ ff ut = p ■ /“"J. Then the solution of the merge 

model is defined to be the projection of R onto 11: that is, 

(/T*. 0 = argmin d [(/“", 0*1 ( 434 > 

(/r'./fben 

where d [(/ 4 " ( , / 5 ol '0> R \ denotes the Euclidean distance between the point (J% ut , f£ ut ) and the ray R: 

d f° 5 u \ R] - mm ||(/ 4 °" r , fn ~ (*. y)|| 2 


5 More generally, as pointed out by|Coclite et al.] j2005| , when the number of incoming links exceeds the number of outgoing links, (Al) and 
(A2) combined are not sufficient to ensure a unique solution. 
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4.2. Well-posedness of the diverge junction model 

In this section, we investigate the well-posedness of the diverge junction model. Unlike previous studies (Garavello 
and Piccoli 2006; Han et al. 2015b[ >, a major challenge in this case is to incorporate drivers’ route choices, expressed 
by the path disaggregation variable p, into the model and the analysis. In effect, we need to establish the continuous 


dependence of the model on the initial/boundary conditions in terms of both p and p. As we shall see in Section 4.2.1 


below, such a continuity does not hold in general. Following this, Section 4.2.2 provides sufficient conditions for the 


continuity to hold. These sufficient conditions are crucial for the desired continuity of the delay operator. 

4.2.1. An example of ill-posedness 

It has been shown by Han et al. ( 2015b] ) that the diverge model (Section 4.1.1 1 with constant turning ratios, a 12 
and a| 3 , is well-posed. However, the assumption of fixed and exogenous turning ratios does not hold in DNL models; 
see ( |3.28| >. As a result the well-posedness is no longer true, which is demonstrated by the following counterexample. 

We consider the diverge junction (Figure |4j and assume the same fundamental diagram /(■) for all the links for 
simplicity. We consider a series of constant initial data parameterized by e on the three links 7), / 2 , and Ip. 


pi(0, x) = pi e ( p c , p jam ), p 2 ( 0 , x) = pf e (p e , p jam ], p 3 ( 0 , x) = pf e (0, p c ) 

where p c and f> jam are the critical density and the jam density, respectively. p\ , pf and pf satisfy: 

/(pf) = sf (pi), /(pf) = (1 - e)/ (pi) 


(4.35) 


(4.36) 


where e > 0 is a parameter. This configuration of initial data implies that link / and link / 2 are both in the congested 
phase, while link 3 is in the uncongested (free-flow) phase, see Figure[6]for an illustration. 


Link /, 



fip) 



Figure 6: An example of ill-posedness of the initial value problem with endogenous route choices. Left: junction topology. Right: illustration of 
the constant initial densities on the links. 


Two paths exist in this example: pi = {/, / 2 ) and p 2 = {/, / 3 ). The initial conditions for p Pl and pf 2 , which arise 
from travelers’ route choices, are as follows. 


p p '(0, x) = e, pf 2 (0, x) = 1-e 


Notice that these initial conditions are part of the initial value problem at the junction. The solution of this initial value 
problem depends on the value of e; in particular, we have the following two cases. 


• When e > 0, we claim that the initial conditions f \, pf and pf satisfying ( |4. 35 [ i-( p4.36] ) constitute a constant 
solution at the junction. To see this, we follow the junction model ( |4.32| ) and the definitions of demand and 
supply \2.6\ - \2.l\ to get 


/r'w 


min 



S 2 U) S 3 (0 ] 

£* 1,2 ’ £* 1,3 / 


= min 



/(Pf) 

E 


c 

1 -e 


£* 1 , 2 -/r« = e /(pt) = /(pf) 

£*1,3-/rw = (l-e)/(Pt) = /(Pf) 


minjc, /(pi). 


C 

1 -e 


- /(Pi) (4.37) 
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(4.39) 
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where C denotes the flow capacity. Thus pi, p® and p| form a constant solution at the junction. 

• When e = 0, the turning ratios satisfy a x ^(t) = 0, o- 1 . 3 (f) = 1. Effectively, link I\ is only connected to link / 3 . 
We easily deduce that 


/f"« = c 

(4.40) 

o 

III 

s: 

(4.41) 

o 

hi 

R 

(4.42) 


As a result, the solution on link I\ is given by a backward-propagating rarefaction wave (or expansion wave, 
fan wave) with p\ and p c on the two sides. On link 1 3 , a forward-propagating rarefaction wave with p c and p? 
(that is, the limit of p^ as e —» 0) on the two sides is created. Link h remains in a completely jam state with full 
density p , “"'. See Figure[7]for an illustration of the solutions. 


I2 




£ > 0 



X 


P2(0, x) = p\ 



0 Plit, 02) S p\ 


t 


X 


b 3 


p 3 ( 0 , x) = p% 


a 3 



0 P3(t, a 3 ) = p% 


t 


£ = 0 



Figure 7: Comparison of solutions on 1 1 , I 2 and I 3 for the two cases: s > 0 and s = 0. The dashed lines represent the kinematic waves 
(characteristics), and darker color indicates higher density. Notice that jumps in the boundary conditions across these two cases exist on I\ and / 3 . 


The two sets of solutions corresponding to s > 0 and = 0 are shown in Figure [7] We first notice that in these 
two cases, the boundary flows are very different, despite the infinitesimal difference in the PDV, namely from s > 0 
to s - 0. In particular, both f° u, (t ) and /™(f) jump from f(p\) and f(pf) (when e > 0) to C (when e = 0). This is a 
clear indication of the discontinuous dependence of the diverge junction model on its initial conditions. 

Let us zoom in on the mechanism that triggers such a discontinuity. According to ( |4.37| ), the expression for 
when e > 0 is 


fi(t) = min C, 


m) c 


i 


min < C, 


g/(Pi) | 


As long as e is positive, the fraction SliEil - f(p x ) < C. However, when e = 0 we have an expression of which 
should be equal to oo since link I 2 has effectively no influence on the junction, and we have min{C, [j) = C. The above 
argument amounts to the following statement: 


if £ > 0 , 
if s = 0 , 


C > 


s/(Pi), 

£ 


C < 


s/(pi) 

£ 


which explains the jump in the solutions when e tends to zero. 
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Remark 4.1. The fact that p e 2 tends to p jam as e —> 0 plays a key role in this example. As we shall see later in Theorem 
14.21 bounding p% away from the jam density is essential for the well-posedness. 


4.2.2. Sufficient conditions for the well-posedness of the diverge model 


As we have previously demonstrated, the diverge model with time-varying vehicle turning ratios may not depend 
continuously on its initial conditions, which are defined in terms of the two-tuple (p, p). In this section, we propose 
additional conditions that guarantee the continuous dependence with respect to the initial data at the diverge junction. 


Our analysis relies on the method of wave-front tracking (Bressan 

(2000 Garavello and Piccoli 2006) and the tech- 

nique of generalized tangent vectors (Bressan 1993; Bressan et al. 

2000 

). In order to be self-contained while keeping 


our presentation concise, we move some general background on these subjects and essential mathematical details to 


Appendix A 

Theorem 4.2. (Well-posedness of the diverge model) Consider the diverge junction (Figure^ and assume that 

1. there exists some 6 > 0 such that the supplies S 2 (t) > 6, S 3 (f) > 6 for all t; 

2. the path disaggregation variable p(t, x) has bounded total variation in t. 

Then the solution of the diverge junction depends continuously on the initial and boundary conditions in terms of p 
and p. 


Proof. The proof is long and technical, and thus will be presented in Appendix B. 1 following necessary mathematical 
preliminaries in Appendix A □ 


4.3. Well-posedness of the merge model 

For the merge junction depicted in Figure |4j the path disaggregation variables p become irrelevant since there is 
only one downstream link. According to Theorem 5.3 of Han et al. ( 2015b) >, the solution at the merge junction, in 
terms of density p, depends continuously on the initial and boundary conditions. Moreover, according to ( |3.9| ), the 
propagation speed of p is the same as the vehicle speed v(p), we thus conclude that p also depends continuously on 
the initial and boundary conditions. This shows the well-posedness of the merge junction. 


5. Continuity of the delay operator 

The proof of the continuity of the path delay operator is outlined as follows. We first show that the first hypothesis 


of Theorem 4.2 holds for networks that consists of only the merge and diverge junctions discussed in Section |4. 1 [ this 
will be done in Section [5TT| Then, we propose some mild conditions on the fundamental diagram and the departure 


rates in Section 5.2 to ensure that the second hypothesis of Theorem 4.2 holds. It then follows that both the merge and 


diverge models are well posed. Finally, in Section 5.3 we show the well-posedness of the model in terms of p and p 
at each origin, where a point queue may be present. Put altogether, these results lead to the desired continuity for the 
delay operator. 


5.1. An estimation of minimum network supply 

This section provides a lower bound on the supply, which is a function of density, on any link in the network 
during the entire time horizon [0, T ]. Our finding is that the jam density can never be reached within any finite time 
horizon T , and the supply anywhere in the network is bounded away from zero. As a consequence, grid lock will 
never occur in the dynamic network. 

We denote by D the set of destinations in the augmented network with virtual links (see Section [T3| : that is, every 
destination d £ D is incident to a virtual link that connects d to the rest of the network; see Figure [ 8 ] For each d e 1). 
we introduce its supply, denoted S' 1 . to be the maximum rate at which cars can be discharged from the virtual link 
connected to d. Effectively, there exists a bottleneck between the virtual link and the destination; and the supply of 
the destination is equal to the flow capacity of this bottleneck; see Figure [ 8 ] Notice that in some literature such a 
bottleneck is completely ignored and the destination is simply treated as a sink with infinite receiving capacity. This 
is of course a special case of ours once we set the supply S d to be infinity. However, such a supply may be finite and 
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Bottleneck 
Virtual link 


(d) Destination 


Figure 8: Illustration of the supply (receiving capacity) of the destination d. 


even quite limited under some circumstances due to, for example, ramp metering, limited parking spaces, or the fact 
that the destination is an aggregated subnetwork that is congested. 

We introduce below a few more concepts and notations. 


L = min L, : the minimum link length in the network, including virtual links. 

iium 

C mm = min C, : the minimum link flow capacity in the network, including virtual links. 
A = max f,'(p jam -)\ : the maximum backward wave speed in the network. 

/jSJ^ I ' 1 


p = min | p 1 , 1 - //} > 0 , where p J is the priority parameter for the merge junction J (see Section 
denotes the set of merge junctions in the network. 


4.1.21, 'Vm 


6 ° = tiling > 0 

deD 


the minimum supply among all the destination nodes. 


6k = min inf inf S , (pft, x)) 

I,eA t<=\k^,(k+l)j\ ^[a u bi\ 


: the minimum supply at any location in the network during time interval 


[kj, ( k + l)q), where k — 0, 1, ... 


Theorem 5.1. Consider a network consisting of only merge and diverge junctions as shown in Figure [7] Given any 
vector of path departure rates, the dynamic network loading procedure described by the PDAE system yields the 
following property in the solution: 


6 k > p k min {<5® , pC mi "} Vfc = 0, 1, ... 


(5.43) 


Proof We postpone the proof until Appendix B.2 for a more compact presentation. 


□ 


Theorem |5 .1 1 guarantees that the supply values anywhere in the network during the period [k^, (k + 1)^] are 
uniformly bounded below by some constant that depends only on k, although such a constant may decay exponentially 
as k increases. Given that any dynamic network loading problem is conceived in a finite time horizon, we immediately 
obtain a lower bound on the supplies, as shown in the corollary below. 


Corollary 5.2. Under the same setup as Theorem |5. /| we have 

min inf Si(pft, x» > p L ^ J {<5® , pC min \ Vf e [0, T] (5.44) 

Ii£jl x£[aj,bj] 1 ' 


and 


min inf inf S iipft, xj) > p^^ min (b® , pC mm \ 

Ii&n xa[a u bi\ta[0,T] 1 ' 


where the operator |_-J rounds its argument to the nearest integer from below. 


(5.45) 
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Proof. Both inequalities are immediately consequences of ( |5.43| >. □ 

Remark 5.3. Corollary |5.2| shows that, for any network consisting of the merge and diverge junctions, the supply 
values are uniformly bounded from below at any point in the spatial-temporal domain. In particular, a jam density 
can never occur. Moreover, such a network is free of complete 


5.2. Estimation regarding the path disaggregation variables 


In this section we establish some properties of the path disaggregation variable p, which serve to validate the 


second hypothesis of Theorem 4.2 


Lemma 5.4. Assume that there exists M > 0 and e > 0 so that the following hold: 

1. For all links I„ the fundamental diagrams f(j are uniformly linear near the zero density; more precisely, f'ip) 
is constant for p e [ 0 , e] for all I , € 5A. 

2. Each f(-) has non-vanishing derivative, i.e. \f((p)\ > e, Ip € [0 ,pf™]. 

3. The departure rates {h p (-),p 6 P) are uniformly bounded, have bounded total variation (TV), and bounded 
away from zero when they are non-zero; i.e. TV(h p ) < M and h p (t ) e {0} U [e, M]for all p € P and almost 
every t e [0, T]. 


Then the path disaggregation variables p(t, x) are either zero, or uniformly bounded away from zero. Moreover, they 
have bounded variation. 


Proof. The proof is moved to | Appendix B.3| □ 

Remark 5.5. Notice that assumptions 1 and 2 of Lemma \5.4\ are satisfied by the Newell-Daganzo (triangular) funda¬ 
mental diagrams, where both the free-flow and congested branches of the FD are linear. Moreover, given an arbitrary 
fundamental diagram, one can always make minimum modifications at p — 0 and p = p c to comply with conditions 1 
and 2. Assumption 3 of Lemma |5.4| is satisfied by any departure rate that is a result of a finite number of cars entering 
the network. And, again, any departure rate can be adjusted to satisfy this condition with minor modifications. 


5.3. Well-posedness of the queuing model at the origin with respect to departure rates 

In this section we discuss the continuous dependence of the queues q s (t ) and the solutions p, and p 1 . with respect 
to the path departure rates h p (t), p eP. The analysis below relies on knowledge of the wave-front tracking algorithm 


and generalized tangent vector, for which an introduction is provided in Appendix A 


Following Herty et al. (2007), we introduce a generalized tangent vector (q s ,f p , ff) for the triplet (q s ,p,,p i ) 
where q s e R is a scalar shift of the queue q s (■), i.e. the shifted queue is q s f) + ?/ v , while the tangent vectors of p, and 
p r ' are defined in the same way as in Appendix A.2 The tangent vector norm of q s is simply its absolute value \q s \. 


Lemma 5.6. Assume that the departure rates {h p f), p € P) are piecewise constant, and let f p be a tangent vector 
defined via shifting the jumps in h p (f. Then the tangent vector (q s ,f' p ,^ p p ) is well defined, and its norm is equal to 
that off p and bounded for all times. 


Proof. The proof is postponed until | Appendix B.4| 


□ 


6 A complete gridlock refers to the situations where a non-zero static solution of the PDAE system exists. Intuitively, it means that a complete 
jam p = p Ja "' is formed somewhere in the network and the static queues do not dissipate in finite time. 
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5.4. Final proof of continuity for the delay operator 

At the end of this paper, we are able to prove the continuity result for the delay operator based on a series of 
preliminary results presented so far. 

Theorem 5.7. (Continuity of the delay operator) Consider a network consisting of merge and diverge junctions 
described in Section [7~7| under the same assumptions stated in Lenima \5.4\ the path delay operator, as a result of the 
dynamic network loading model ( |3.21[ >-( |3.31| l, is continuous. 

Proof. We have shown that at each node (origin, diverge node, or merge node), the solution depends continuously on 
the initial and boundary values. In addition, between any two distinct nodes, the propagation speeds of either p-waves 
or p- waves are uniformly bounded. Thus such well-posedness continues to hold on the network level. Consequently, 
the vehicle travel speed v,(p,) for any /, depends also continuously on the departure rates. We thus conclude that the 
path travel times depend continuously on the departure rates. □ 


The assumption that the network consists of only merge and diverge nodes is not restrictive since junctions with 
general topology can be decomposed into a set of elementary junctions of the merge and diverge type ( Daganzoj 
|1995| . In addition, junctions that are also origins/destinations can be treated in a similar way by introducing virtual 
links. 


Here, we would like to comment on the assumptions made in the proof of the continuity. In a recent paper (Bres- 


san; 


and Yu, 2015|), a counterexample of uniqueness and continuous dependence of solutions is provided under certain 


conditions. More precisely, the authors construct the counterexample by assuming that the path disa ggreg ation vari¬ 
ables p 1 ’ have infinite total variation. This shows the necessity of the second assumption of Theorem 1 


4.2 


Moreover, 

again in Bressan and Yu (2015), the authors prove that for density oscillating near zero the solution may not be unique 


(even in the case of bounded total variation), which shows the necessity of the third assumption in Lemma 5.4 


Remark 5.8. Szeto ( 20031 provides an example of discontinuous dependence of the path travel times on the path 
departure rates using the cell transmission model representation of a signal-controlled netw’ork. In particular, the 
author showed that when a queue generated by the red signal spills back into the upstream junction, the experienced 
path travel time jumps from one value to another. This, however, does not contradict our result presented here for 


the following reason: the jam density caused by the red signal in Szeto\ {2003) does not exist in our network, which 
has only merge and diverge junctions (without any signal controls). Indeed, as shown in Corollary \5.2\ the supply 
functions at any location in the network are uniformly bounded below by a positive constant, and thus the jam density 
never occurs in a finite time horizon. 

The reader is reminded of the example presented in Section ed. 2.1\ where the ill-posedness of the diverge model is 
caused precisely by the presence of a jam density. The counterexample from \Sze to\ (2003\ is constructed essentially in 
the same way as our example, by using signal controls that create the jam density. 


We offer some further insights here on the extension of the continuity result to second-order traffic flow models. 
To the best of our knowledge, second-order models on traffic networks are mainly based on the Aw-Rascle-Zhang 
model (Garavello and Piccoli 2006 Haut and Bastin 2007] Herty and Rascle 2006| > and the phase-transition model 
(Colombo et al. 2010[ ). Solutions of the Aw-Rascle-Zhang system may present the vacuum state, which, in general, 
may prevent continuous dependence ( |Godvik and Hanche-Qlsen||2008| l. On the other hand, the phase transition model 
on networks behaves in a way similar to the LWR scalar conservation law model. However, continuous dependence 
may be violated for density close to maximal; see Colombo et al. ( 2010| >. Therefore, extensions of the continuity 
result to second-order models appear possible only for the phase-transition type models with appropriate assumptions 
on the maximal density achievable on the network, but this is beyond the scope of this paper. 


6. Conclusion 


This paper presents, for the first time, a rigorous continuity result for the path delay operator based on the LWR 
network model with spillback explicitly captured. This continuity result is crucial to many dynamic traffic assignment 
models in terms of solution existence and computation. Similar continuity results have been established in a number of 
studies, all of which are concerned with non-physical queue models. As we show in Section 4.2.1 the well-posedness 
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of a diverge model may not hold when spillback occurs. This observation, along with others made in previous studies 
(Szeto 2003| , have been the major source of difficulty in showing continuity of the delay operator. In this paper, we 
bridge this gap through rigorous mathematical analysis involving the wave front tracking method and the generalized 
tangent vectors. In particular, by virtue of the finite propagation speeds of p-waves and p-waves, the continuity of 
the delay operator boils down to the well-posedness of nodal models, including models for the origins, diverge nodes, 
and merge nodes. Minor assumptions are made on the fundamental diagram and the path departure rates in order 
to provide an upper bound on the total variations of the density p and the path disaggregation variables p, which 
subsequently leads to the desired well-posedness of the nodal models and eventually the continuity of the operator. 

A crucial step of the above process is to estimate and bound from below the minimum network supply, which is 
defined in terms of local vehicle densities. In fact, if the supply of some link tends to zero, the well-posedness of 
the diverge junction may fail as we demonstrate in Section 4.2.1 This has also been confirmed by an earlier study 


(Szeto 20031, where a wave of jam density is triggered by a signal red light and causes spillback at the upstream 
junction, leading to a jump in the path travel times. Remarkably, in this paper we are able to show that (1) if the 
supply is bounded away from zero, then the diverge junction model is well posed; and (2) the desired boundedness 
of the supply is a natural consequence of the dynamic network loading procedure that involves only the simple merge 
and diverge junction models. This is a highly non trivial result because it not only plays a role in the continuity proof, 
but also implies that gridlock can never occur in the network loading procedure in a finite time horizon. However, we 
note that in numerical computations gridlock may very well occur due to finite approximations and numerical errors, 
while our no-gridlock result is conceived in a continuous-time and analytical (i.e., non-numerical) framework. 

It is true that, despite the correctness of our result regarding the lower bound on the supply, zero supply (or 
gridlock) can indeed occur in real-life traffic networks as a result of signal controls. In fact, this is how Szeto (2003 


constructs the counter example of discontinuity. Nevertheless, signal control is an undesirable feature of the dynamic 
network loading model for far more obvious reasons: the on-and-off signal control creates a lot of jump continuities 
in the travel time functions, making the existence of dynamic user equilibria (DUE) almost impossible. Given that the 
main purpose of this paper is to address fundamental modeling problems pertinent to DUEs, it is quite reasonable for 
us to avoid junction models that explicitly involve the on-and-off signal controls. Nevertheless, there exist a number 
of ways in which the on-and-off signal control can be approximated using continuum (homogenization) approaches. 
Some examples are given in Aboudolas et al.| (|2009)>; |Gazis and Potts|(|1963|);|Han et al.|(|2014|l and|Han and Gay ah 
J20I5) . 

Regarding the exponential decay of the minimum supply, we point out that, as long as the jam density (gridlock) 
does not occur, our result is still applicable. From a practical point of view most networks reach congestion during 
day time but are almost empty during night time. In other words, the asymptotic gridlock condition is never reached 
in reality. Therefore, we think that our result is a reasonable approximation of reality for DUE problems. 


We would also like to comment on the continuity result for other types of junction models. Garavello and Piccoli 


( |2Q06j ) prove that Lipschitz continuous dependence on the initial conditions may fail in the presence of junctions with 
at least two incoming and two outgoing roads. Notably, the counterexample that they use to disprove the Lipschitz 
continuity is valid even when the turning ratio matrix A J is constant and not dependent on the path disaggregation 
variables It is interesting to note that, in the same book, continuity results for general networks are provided, but 
only for Riemann Solvers at junctions that allow re-directing traveling entities - a feature not allowed in the dynamic 
network loading model but instead is intended for the modeling of data networks. 

The findings made in this paper has the following important impacts on DTA modeling and computation. (1) The 
analytical framework proposed for proving or disproving the well-posedness of junction models and the continuity 
property can be applied to other junction types and Riemann Solvers, leading to a set of continuity/discontinuity results 
of the delay operators for a variety of networks. (2) The established continuity result not only guarantees the existence 
and computation of continuous-time DUEs based on the LWR model, but also sheds light on discrete-time DUE 
models and provides useful insights on the numerical computations of DUEs. In particular, the various findings made 
in this paper provide valuable guidance on numerical procedures for the DNL problem based on the cell transmission 
model (Daganzo 1994 1995), the link transmission model (Yperman et al. 2005), or the kinematic wave model (Han 


et al.| 2015b Lu et al. 2013)>, such that the resulting delay operator is continuous on finite-dimensional spaces. As 


result, the existence and computation of discrete-time DUEs will also benefit from this research. 
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Appendix A. Essential mathematical background 

Appendix A.l. Wave-front tracking method 

The wave-front tracking (WFT) method was originally proposed by Dafermos ( 1972[ > as an approximation scheme 
for the following initial value problem 


J dtp(t, x) + d x f(p(t, x)) = 0 
}p( 0 , x) = p(x) 


(A.l) 


where the initial condition p(-) is assumed to have bounded variation (BV) (Bressan 2000| ). The WFT approximates 
the initial condition p(-) using piecewise constant (PWC) functions, and approximates /(•) using piecewise affine 
(PWA) functions. It is an event-based algorithm that resolves a series of wave interactions, each expressed as a 
Riemann Problem (RP). The WFT method is primarily used for showing existence of weak solutions of conservation 
laws by successive refinement of the initial data and the fundamental diagram; see Bressan ( 2000) > and Holden and 
Risebro (1995 ). |Garavello and Piccoli| ( |2006) ) extend the WFT to treat the network case and show the existence of the 

weak solution on a network. We provide a brief description of this procedure below. 

Fix a Riemann Solver (RS) for each road junction in the network. Consider a family of piecewise constant 
approximations p' ; (x) of the initial condition p,(x) on each link and a family of piecewise affine approximation of the 
fundamental diagrams ffi(p), where s is a parameter such that p?(x) —> p,(x) and ffi(p) —* f(p) as s —> 0. A WFT 
approximate solution on the network is constructed as follows. 

1. Within each link, solve a Riemann Problem at each discontinuity of the PWC initial data. At each junction, 
solve a Riemann Problem with the given RS. 

2. Construct the solution by gluing together the solutions of individual RPs up to the first time when two traveling 
waves interact, or when a wave interacts with a junction. 

3. For each interaction, solve a new RP and prolong the solution up to the next time of any interaction. 

4. Repeat the processes 2-3. 

To show that the procedure described above indeed produces a well-defined approximate solution on the network, one 
needs to ensure that the following three quantities are bounded: ( 1 ) the total number of waves; ( 2 ) the total number of 
interactions (including wave-wave and wave-junction interactions); and (3) the total variation (TV) of the piecewise 
constant solution at any point in time. These quantities are known to be bounded in the single conservation law case; 


in fact, they all decrease in time (Bressan 2000). However, for the network case, one needs to proceed carefully in 


estimating these quantities as they may increase as a result of a wave interacting with a junction, which may produce 
new waves in all other links incident to the same junction. The reader is referred to Garavello and Piccoli ( j2006[ ) for 
more elaborated discussion on these interactions. For a sequence of approximate WFT solutions p e , s > 0, if one can 
show that the total variation of p e is uniformly bounded, then as s —> 0 a weak entropy solution on the network is 
obtained. 

Appendix A.l. Generalized tangent vector 

The generalized tangent vector is a technique proposed by Bressan ( 1993} ) to show the well-posedness of con¬ 
servation laws, and is used later by |Garavello and Piccoli (2006 1 to show the well-posedness of junction models in 
connection with the LWR model. Its mathematical contents are briefly recapped here. 

Given a piecewise constant function F(x) : [a, b] —> R, a tangent vector is defined in terms of the shifts of the 
discontinuities of F(-). More precisely, let us indicate by {vt}^ the discontinuities of F where a - xq < X\ < ■ ■ ■ < 
xn < xjv+i = b , and by {F^f * 1 2 3 4 the values of F on (jc*_i, xf). A tangent vector of F(f) is a vector g = (gi,..., g^) R. N 
such that for each e > 0, one may define the corresponding perturbation of F(f), denoted by /-’' (•) and given by 

F e (x) = F k x e [x k -\ + ef k -i, x k + eg k ) 

for k — 1,.... N + 1, where we set go = ftv+i = 0. The norm of the tangent vector g is defined as 


ll£ll = Y^\fk\-\F{x k +)-F{x k -)\ 


(A.2) 


k= 1 
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In other words, the norm of the tangent vector is the sum of each \£k\ multiplied by the magnitude of the shifted jumps. 

The procedure of showing the well-posedness of a junction model using the tangent vectors is as follows. Given 
piecewise constant initial/boundary conditions on each link incident to the junction, one considers their tangent vec¬ 
tors. By showing that the norms of their tangent vectors are uniformly bounded in time among all approximate 
wave-front tracking solutions, it is guaranteed that the L 1 distance of any two solutions is bounded, up to a multiplica¬ 
tive constant, by the L l distance of their respective initial/boundary conditions. More precisely, we have the following 
theorem 


Theorem Appendix A.l. If the norm of a tangent vector at any time is bounded by the product of its initial norm 
and a positive constant, among all approximate wave-front tracking solutions, then the junction model is well posed 
and has a unique solution. 

Proof See Garavello and Piccoli ( f2006| l. P 


Throughout this paper, we employ the notation (p ; , pj) to represent a wave interacting with the junction from road 
/,, where pj is the density value in front of the wave (in the same direction as the traveling wave) and p, is the density 
behind the wave. After the interaction, a new wave may be created on some road f (it is possible that f = If), and we 
use pj and pj to denote the density at J on road 7 ; before and after the interaction, respectively. 


4.1.1 


Lemma Appendix A.2. Consider the diverge model in Section 
the shift produced on any Ij, as a result of the shift ft on /„ satisfies 


If a wave (p;, p ( - ) on f interacts with J then 


where 

A<?, - fi(pt)-fi(pj), 

are the jumps inflow across these two waves. 


A g; 

A <y, 


fdpj 


A qj = fj(pj)-fj(pf) 


(A3) 


Proof. To fix the ideas, we assume that /, is the incoming road f (the other cases can be treated in a similar way). Let 
(p\,pj) be the interacting wave, then it must be that p\ < cr. 

If no wave is reflected on road f after the interaction (i.e. no new backward wave (p \, p j) is created on f as 


a result of the interaction), then we can conclude the lemma by following the same estimate as shown in Garavello 
and Piccoli (2006 1 . If, on the other hand, a wave (pi,p|) is reflected then pj < p\, (p\,pj) is a rarefaction wave. 


and f(pj) < f(pj) < pi ■ In other words, part of the change in the flow caused by the interaction passes through the 
junction, and part of it is reflected back onto I\. We may then apply the same argument as in Garavello and Piccoli 
(2006) to the part that passes through the junction (f(p \) - f(pj)) and conclude the lemma. □ 


Definition Appendix A.3. If a p-wave from /, interacts with the junction, producing a p-wave on link Ij, we call I , 
and Ij the source and recipient of this interaction, respectively. Such an event is denoted f —* I j. 


We observe that \£i(pj - pi)\ is precisely the L 1 -distance of two initial conditions on /, with and without the shift f; 
at the discontinuity (p,-, pj). Similarly, | fjipj - pj)\ is the L 1 -distance of the two solutions on Ij as a result of having 
or not having the initial shift £; on /,-, respectively. Furthermore, if is the only shift in the initial condition on then 
| £i(pj - pj) | is the norm of the tangent vector for /, (see definition < |A .2[ )) before the interaction, and fijpj - pj)\ is the 
norm of the tangent vector for Ij after the interaction. From (|A.3[), we have 


\£M -Pi)' 


\Aqj\ 

|A?/| 


1 Z(p7 


Pi)| 


(A.4) 


Therefore, to bound the norm of the tangent vectors one has to check that the multiplication factors j^j remain 
uniformly bounded, regardless of the number of interactions that may occur at this junction. 

Let us consider the diverge junction, and assume that a wave interacts with J from road /, and produces a wave 
on road Ij, where i, j — 1,2,3. As in the proof of Lemma Appendix A.2 we can restrict to the flow passing through 
the junction (indeed the reflected part has the same flow variation and smaller shift since the reflected wave is a slow 
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big shock). According to ( | 4 . 32 | >, Ar/ 2 = 042 Agi and A<73 = <343 A/4 . Consequently, we have the following matrix of 
multiplication factors j^j: 


|Ar yi | 

|Ar/i| 

= 1, 

IA92I 

IA91I 

- 0'1,2» 

1A93I 

ia?ii 

- a 1,3; 

|A< 7,1 

1 

IA92I 

= 1 , 

1A931 

0:1,3 

IA92I 

a 1,2’ 

IA92I 

|A? 2 | 

a 1,2 

|A< 7,1 

1 

IA92I 

a i f 2 

IA93I 

= 1. 

IA93I 

<*1,3 ’ 

IA93I 

<*1,3 ’ 

IA93I 


(A.5) 


We denote this matrix by {<2i/};j=i,2,3- According 1 
vector norm, it suffices to keep tract of just one single shift and show the corresponding tangent vector norm is 
bounded regardless of the number of wave interactions. To this end, we consider the only meaningful sequence of 

wave interaction, which is of the form /, —> / ; , /, —> 4,_We observe that, if a 1,2 and a 13 are nonzero constants, 

then QqQjk = Qik for any i, j, k = 1,2,3. This means that no matter how many interactions occur, the multiplication 
factor is always an element of the matrix {Qij} and is uniformly bounded. Thus the diverge model with fixed turning 
ratios is well posed. However, in the DNL model where a\i and a 13 are time-varying, this is no longer true as we 


have shown in the counterexample in Section 4.2.1 The well-posedness requires some additional conditions to hold. 


and the corresponding proof needs more elaborated arguments that take into account the p- waves. These will be done 
in Theorem l4.2l 


Appendix B. Technical Proofs 

Appendix B.l. Proof of Theorem \ 4. 2\ 

Proof. This proof is completed in several steps. 


[Step 1]. We invoke the wave-front tracking framework and the generalized tangent vector to show the desired con¬ 
tinuous dependence. As discussed in Appendix A.2 it boils down to showing that the increase in the tangent vector 
norm, as a result of arbitrary number of wave interactions (including both p -wave and p- wave), remains bounded. 


[Step 2]. First of all, the end of |Appendix A. 2] shows that, for constant values of p (that is, with fixed vehicle 
turning ratios), the tangent vector norm is uniformly bounded regardless of the interactions between the /(-waves and 
the junction. Next, one needs to consider the case where p changes from one value to another, i.e. when a p -wave 
interacts with the junction. In general, we consider two consecutive interaction times of the p- waves: 4 and 4+1, and 
let /3k be the multiplication factor for the tangent vector norm of p relative to the times 4 and 4+1. More precisely, we 
have v(4+i~) = /4-v(4+X where v(t) is the tangent vector norm of p at time f. Clearly, /4 can only take value in the 
following matrix where a\ ^ and or*, are given by the constant p value during (4, 4+1); see ra- 


Q = {Q-j) = 


1 

a \,2 

(Y^ 

a \,3 

l/< 2 

1 

a\ z /a\ 

i/< 3 

a i, 2 /°U ,3 

1 


k = 0, 1, 2, 


(B.l) 


Similar arguments can be carried out to all other intervals [4+1, 4+2], [4+2, 4+3] and so forth. In the rest of the 
proof, our goal is to estimate the multiplication factor for the tangent norm of p under the WFT framework, assuming 
arbitrary interaction patterns of the /(-waves and p- waves. 


[Step 3]. We define y = 4- > 0 where <5 is stated in the hypothesis of this theorem and C\ is the flow capacity of link 
1 

l\. Without loss of generality, we let y < Throughout Step 3, we assume that 0-12(4 > y and 043(f) > y (the other 
cases will be treated in Step 4). 
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As demonstrated at the end of Appendix A.2 in order to estimate the increase in the tangent vector norms it 
suffices to consider the following type of sequence in which p -waves interacts with the junction: /, —> /, then I j —> 4; 
in other words, the recipient of the previous p -wave interaction is the source of the next p-wave interaction. For each 
k > 1, we let I ® e {/[, 1 2 , 4) be the recipient of the last interacting p-wave during (4, 4+i). This gives rise to the 
sequence see Figure | bT 9] for an illustration. 


x p -wave interaction 
O Last p -wave interaction 


7 -(*) 


r(*+i) 


r(.k+ 2) 


-X-X- 


O • 


l k+l 


l k+2 


' h+3 


O • x —x—x—x— x o • 




o k a k 

'-^'1 2 5 ^1 3 


a k +1 a k+1 
m 1,2 ’^ 1,3 


a k+2 sy k+ 2 


Figure B.9: t^s represent the times when //-waves interact with the junction, changing the turning ratios ai t 2 and a\ t 3 . The crosses represent the 
times at which p-waves interact with the junction. The circles indicate the last p-wave interaction within the interval during which p is constant. 
/W’s denote the recipients of the last interacting p-waves. 


We make the following crucial observation. Consider any three elements in the sequence of the form I ik) = /,, 
f (k+ 1 ) = 4 and I (k+2 > = Jj where i, j e {1, 2, 3). By definition, the product of the multiplication factors within 
(4+t, 4+3) is 


Qa l ' Q\f - A\ +l ■ Q‘j where A k+] e 


(yk+l 

a l,2 

ry k+ 1 
a l,2 

a k+2 

a l,3 

(Y k+ 2 
a \,2 

a k+2 

a \,2 

ry k+ 1 
^1,3 


le{k+l , k + 2] 


(B.2) 


where we use the superscripts to indicate dependence on a specific time interval. The significance of ( |B.2| ) is that the 
multiplication factor can be decomposed into a term A 2 with a very specific structure (to be further elaborated below), 
and a term that would have been the multiplication factor as if the middle link / 2 = l <k+l> were removed. Clearly, the 
argument above applies equally to the case where I tk * 11 = 4, and we use 


A 


r-+i 

3 



ry k+ 1 
^1,3 

a k+2 

a l,3 


yk +1 

J +2 
l \,3 


a v. 2 2 ) 

a U2 } 


to represent the term factored out if the middle link I ' k+1 ’ = 4 is removed. By repeating this procedure, one may 
eliminate all links 4 and 4 from the sequence {/®} except when they are the first or possibly the last in this sequence. 
As a result, the multiplicative terms (A| +1 } and |,4* 11 ) are factored out. Therefore, the entire multiplication factor for 
the tangent norm is bounded from above by 


FK*' U A 


k „+1 
3 


m= 1 


n= 1 


(B.3) 


where {4n)” =1 and {k„}“ j are subsequences of { k } such that 4 im+1) = 4 and /L+D = / 3 xhe first two terms in ( |B.3| > 
are results of eliminating 4 and 4 from the sequence {l (k> \; the third term takes care of the first element and possibly 
the last element in the sequence, and is derived from the observation that (4 < i for all i, j e {1, 2, 3} and k > 0. 

It remains to show that the first and second terms of ( |B.3| > are bounded, and hence the entire ( |B.3[ ) is bounded. We 
write, without referring explicitly to the multiplication indices, that 


OO 


n+L'. 

m= 1 



k+2 

a l,3 


a 


k +1 
1,3 
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This leads to the following. 


fl4'” +1 = ex P X loga\i ~ l°S a u 2 ex P ^og(l - a\f) - log(l - a^) 


k 


exp 


SU P , , 

ae[y, l-y] \ a \ 

1 




V k 
\ 


*+1 fyk+2\ 

a \,2 I 


• exp 


sup 

^are[y, l-y] 11 — 




k+2 tr k+l I 

“1,2 I 


r 


< exp - ■ TV(i-i) ■ exp - ■ TV(p) < oo 


1 


r 


where TV(p) is the bounded total variation of the path disaggregation variable p. 


[Step 4]. In this part of the proof we deal with the situation where ct\g(t) < y\ the other case where a\p(t) < y is 
entirely similar. We begin with the following observation. 


/r« = min 



S 2 {t) S 3 (f) 1 
O'!.2W’ ari i3 (f)J 


min 


Diit) 


S 3 jt) \ 
a\,iit) ] 


since = s 2 (t>c I > q > Dft). In other words, the minimum is never attained at S2<n . and any p- 

wave interacting from P will not generate new waves on I\ or / 3 since does not change before and after the 

interaction. So the only interactions that may change the tangent vector norm are 7/ —> / 3 and / 3 —» I\, and the 
corresponding multiplication factors are respectively ai j3 (0 and 

Consider any sequence of interaction times {4} of the /j-waves. Since the interactions I\ —7 /] or / 3 —> / 3 do 
not generate any increase in the tangent norm, we only need to consider the following sequence of interactions: 
... —> I\ —> / 3 —> / [ —> / 3 _The resulting multiplication factor is bounded by constant of the form 



where the second multiplicative term deals with the first and the last element in the sequence since max{aj 3 , -!-} = 
-j- < pL. f or a ]j £ We may then proceed in the same way as Step 3 and conclude that this multiplication factor is 
bounded, provided that TV(p) <00. □ 


Appendix B.2. Proof of Theorem \5. 1\ 

Proof. Since 5k is defined in terms of the supplies, it suffices for us to focus only on those densities such that 
pft, x) > p' for f 6 3K, ( t, x) e [0, T] x [a,-, /?,]. It is also useful to keep in mind that densities beyond the criti¬ 
cal density always propagate backwards in space. The proof is divided into several steps. 

[Step 1], k = 0. We have t e [0, j). Since the network is initially empty, all the supplies S fpf t, x)), /,■ 6 (f, x) 6 

{0)x[a,. bi], are maximal and equal to the respective flow capacities at t = 0. Afterwards, a higher-than-critical density 
or a lower-than-maximum supply can only emerge from the downstream end of a link and propagate backwards along 
this link. Moreover, these backward waves can never reach the upstream end of the link within [0, ) since is the 

minimum link traversal time for backward waves. 

A higher-than-critical density (backward wave) can only arise in one of the following cases. 

Case (1). A forward wave from I\ (see Figure Q interacts with the diverge junction and creates a backward wave on 

h- 

Case (2). Aforw’ard wave from I 4 (or f) interacts with the merge junction and creates a backward wave on I 4 or f. 
Case (3). Aforw’ard wave interacts with a destination from the relevant virtual link and creates a backward wave on 
the same virtual link. 

Each individual case will be investigated in detail below. 
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• Case (1). We use the same notation shown in Figure [4] According to the reason provided above, S ,(t) = C,, 
i = 2,3 for t e [0, k). Let the time of the wave interaction be t. Recall from (|4.32|) that 


= min 



_Ci_ C 3 \ 

<*1,2© ’ at,3©j 


(B.4) 


If the minimum is attained at D\(i), then the entrance of 7) will remain in the uncongested phase, i.e. p \(/+, b\—)< 
Pj after the interaction. Hence, no lower-than-maximum supply is generated. On the other hand, if say j s 
the smallest among the three, then 

^(P!(/+,/?!-)) > /i(p,(?“+, &,-)) = f{ ,ut (t+) = —%- > C 2 > C min (B.5) 

at ,i(t+) 


In summary, the supply values corresponding to the backward waves generated at l \, if any, are uniformly above 

£min 

• Case (2). We turn to the merge junction shown in Figure^ and note S f ,(t) = Q, for t g [0, j). As usual, we let 
t be the time of interaction. There are two further cases for the merge junction model as shown in Figure [5] We 
first consider the situation illustrated in Figure |5j a). Clearly, we have that 

s 4 (p 4 (t+, b 4 -)) > f 4 (p 4 (t+, b 4 -)) = f° A u, (t+) = P C 6 > pc min (B.6) 

s 5 (P5(t+, b 5 -y) > f 5 (P5(t+, b s -yy = /r'(F+) - (1 -P)C 6 > pc min (B.7) 


That is, the supply values of the backward waves generated at the downstream ends of I 4 and Is are uniformly 
above pC' mn . 

For the situation depicted in Figure Bib), we first note that the coordinates of Q' is ( pC 6, (1 - p)C(>), and the 
coordinates of the solution Q = (/J“ ’ , either satisfy 

/;“'•* < P C 6 and f'2’* > (1 -p)C 6 (B.8) 

as shown exactly in Figure [5] b), or 

C’* > pC 6 and / 5 ou '-* < Cl - p)C 6 (B.9) 

Taking ( |B.8[ > as an example (the ( ]B.9[ i case can be treated similarly), we have 

s 5 (Ps(t+, bs-yy > fs(p$(.t+, bs-yy > (i -p)c 6 > pc min ■ cb.io) 


and no increase in density is produced at the downstream end of link I 4 since its exit flow is equal to the demand. 

To sum up, the supply values corresponding to the backward waves generated at I 4 or Is, if any, are uniformly 
above pC min . 

• Case (3). For each destination d e D, let S d be its supply and denote by vl the virtual link connected to d. If 
S d > C v i then no backward waves can be generated on this link and the supply value on the link is always equal 
to C v /. If S ' 1 < C v i, then only one higher-than-critical density can exist on this virtual link, that is, p such that 
P>Pw and f,i(p ) = S d . Clearly, its supply value is equal to S' 1 . 

To sum up, the supply value corresponding to the backward waves generated at any virtual link, if any, is 
bounded below by 6 °. 


Finally, we notice that all the backward waves exhaustively described above originate from the downstream end 
of a link, and they cannot reach the upstream end of the same link within period [0, j). Thus, these backward waves 
cannot bring further reduction in the supply values by means of interacting with the junctions. We conclude that 
during period [0, j), the minimum supply in the network, do, is bounded below; that is, 

do > min Jd® , pC min J 
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[Step 2]. We move on to k > 1. In addition to Case (1 )-(3), which do not bring any supply values below min[d®, pC mm }, 
two more cases may arise in which higher-than-critical densities may be generated as a result of backward waves in¬ 
teracting with junctions: 

Case (4). A backward wave from h (or f, see Figure^ interacts with the diverge junction and creates a backward 
wave in I\. 

Case (5). A backward wave from f (see Figure^ interacts with the merge junction and creates a backward wave in 
I 4 or f. 

Case (4) and Case (5) are analyzed in detail below. 

• Case (4). Without loss of generality, we assume that the backward wave that interacts with the diverge junction 
is coming from 1 2 , and has the density value pj e (pf pf"‘]. Let t he the time of the interaction. In view of 
dB4l >, if the minimum is attained at D i(f), then the interaction does not bring any increase in density at the 
downstream end of /[, hence no decrease in the supply there. 

If the minimum is attained at we deduce in a similar way as that 

S i(pi(t+, b\-)) > fi(pi(t+, b\-j) = fr\t+) = > S 2 (p 2 ) > tft-i 

The last inequality is due to the fact that a backward wave such as p 1 must be created at the downstream end 
of I 2 at a time earlier than kj, thus its supply value S 2 (pj) must be bounded below by S k 1 . If the minimum is 
attained at ,V ' K 2 , we have 

ai,3© 


S i(p\(t+, bi~)j > /i(pi(f+, bi-j) = f° u \i+) = —> S 3 (f+) > S k - 1 

£*1,30+) 

The last inequality is because any lower-than-maximum supply on h, must be created in the previous time 
interval. 

To sum up, the supply values corresponding to the backward waves generated at l \, if any, are bounded below 
by S k -i. 

• Case (5). For the merge junction, we begin with the case illustrated in Figure [5ja). Assuming the backward 
wave that interacts with the merge junction from / f , has the density value p 6 e (p c 6 , p J ( ‘" n ] ■ Similar to ( |B.6[ )-( p3T7| i, 
we have 


S 4 (p4(t+, b 4 -)) > f 4 (p 4 (t+, b 4 -)) = f 4 u, (t+ ) = pS 6 (p 6 ) > p5 k - 1 (B.12) 

S 5 (P5(J+, b 5 -)) > f 5 (p 5 (t+ , b 5 -)) = f° u '(t+) = (1 -p)S 6 (p 6 ) > p5 k CB. 13) 

where the last inequalities are due to the same reason provided in Case (4). The case illustrated in Figure[5}b) 

can be treated in the same way as ( |B .8[ >-( [BTT01 >. 

To sum up, the supply values corresponding to the backward waves generated at I 4 or f, if any, must be bounded 
below by p5 k -1 ■ 

So far, we have shown that for t e \kk, (k + 1 j-f where k > 1, the presence of higher-than-critical densities, as ex¬ 
haustively illustrated through Case (1)- Case (5), brings supply values throughout the entire network that are bounded 
below by 

minjmin {6°, pC m "’} , S k -\ > P&k-\ 1 = min ^6° , pC m,n , p6 k -\\ = min{(5 s , (B-14) 

Case (l)-(3) Case (4) Case(5) 

[Step 3]. Recall ( |BTTT| ) and ( |BTT4] ): 

do ^ min {d® , pC m '"}, 5 k > minjd® , ^d^ 1 } Vk > 1 (B. 15) 
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We define 6q = min j 6° , pC mm } and 6k = min {b® , p6k- i }, k > 1. Clearly, p6u- i < 6° for all k > 1, which implies 
that 6k = min j<5® , p6k- 1 } = p6k- 1 , VA: > 1. Thus 

6 a > §o = minjd®, pC m "'J, 6k > 6k = p6k- \ = p k 6o - /^minjd® , pC mm } VA: > 1 

n 


Appendix B.3. Proof of Lemma \5.4\ 

Proof The proof is divided into a few steps. 

[Step 1]. Notice that assumption 3 implies that the departure rates {h p (-), p e V) are non-zero on a finite set of inter¬ 
vals. Indeed, let n be the number of intervals where h p {-) does not vanish then ne < TV(h p ) < oo, which implies that 
n is finite. 

[Step 2]. Using the assumption that each fundamental diagram /■(■) has nonvanishing derivative, we have, at each 
origin, that the bounded total variation of the density must imply bounded total variation of the flow, and vise versa. 
Therefore, the assumption on the bounded variation of h p {-) implies that the density and path disaggregation variable 
are of bounded variation on the virtual link incident to that origin. 

[Step 3]. Assumption 1 guarantees that, on each link, all waves of the type (0, p) or ( p , 0) are contact discontinuities 
for p sufficiently small and, in particular, travel with a constant speed. Moreover, all p -waves travel with the same 
speed for low densities. Therefore, taking into account Step 1, whenever p is non-zero, it is uniformly bounded away 
from zero on the entire network. In particular there exists e' > 0 so that p p {t, x) e (0) U [e', 1] for every t, x, i and p. 
Therefore, at diverging junctions, the coefficients cri^O and ac\p,(t) satisfy the same properties. 

[Step 4]. Let us now turn to the total variation of the path disaggregation variables. We know, from Step 2, that p 
has bounded variation on virtual links incident to the origins. We also know from Step 3 that p is bounded away from 
zero whenever it is non-zero; and the same holds for the turning ratios at diverge junctions. Consider a p- wave (p/,p r ), 
then its variation \pi - p,\ can change only upon interaction with diverge junctions. More precisely, denote by (pj,p~) 
and (pfp*) the wave respectively before and after the interaction, we have \pf - p+\ < £\pj - p~\ where a = a\p_ or 
a = a 1,3 and is bounded away from zero. Since the //-waves travel only forward on the links with uniformly bounded 
speed, we must have that the interactions with diverge junctions can occur only finite number of times. Thus we 
conclude that p has bounded total variation. O 

Appendix B.4. Proof of Lemma |5.6| 

Proof Let f p be the shift of the j-th jump of h p (-), which occurs at time f ; , then the expression of q s (f) is possibly 
affected on the interval [tj, tj + C p \ (assuming f p > 0). More precisely, if q s (t) > 0 then no wave is generated on the 
virtual link while a shift in the queue is generated with // v = A jh p ■ f p , where A fh p is the jump in h p (-) occurring at time 
tj. On the other hand, if qf t) = 0 then a wave ip ,p 4 j is produced at time t } on the virtual link with shift £f = A • f J p 
where A is the speed of the wave (p ,(>' ). We can then compute: 

t P ■ (P + ~P~) = & ■ * ■ (P + ~pn = & ■ f(p+) ~ f( _ p - } ■ ( p + - p -) 

1 y p + - p 

Moreover, f(p + ) - f(p ) = Ajh p , thus the norm of the tangent vector generated is the same as before. 

To prove that the norm of the tangent vector (r] s , tf p ) is bounded for all times, let us first consider a backward- 
propagating wave (p~,p + ) interacting with the queue q s , and with a shift f p . Then we can write: 

<t s -<fs= f(P~) ~ .f(P + ) 
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where q s and q + s are the time-derivatives of q s (-) before and after the interaction, respectively. Therefore, denoting A 

the speed of the wave (p~,p + ), we get: 

Vs = %■ (ZOO - /(P + )) = % _ ( f r P ,\ ' (J(P~) ~ f(P + )) = t'p ■ (P~ ~ P + ) 

A f(P )~ HP ) 

Thus the norm of the tangent vector is bounded. 

The norm of the tangent vector (q s ,$‘^ p ) may also change when the queue q s becomes empty, but this can be 
treated in the same way as in |Herty et al.| < (2007] ). □ 
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